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Abstract 

A dynamical model of confinement based on a transport theoretical de- 
scription of the Friedberg-Lee model is extended to explicit color degrees of 
freedom. The string tension is reproduced by an adiabatic string formation 
from the nucleon ground state. Color isovector oscillation modes of a qq-system 
are investigated for a wide range of relative gg-momenta and the dynamical 
impact of color confinement on the quark motion is shown. 
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Kinetic theories are a well established tool in describing intermediate energy 
heavy-ion collisions. Their success has stimulated the development of transport the- 
ories based on QCD to describe collisions for even higher energies. For relativistic 
energies, string and parton cascade models (Fritjof |l|], Venus RQMD HIJING 
1^, parton cascade have been developed. A serious problem of these models is 
the description of the final state of the reaction, the hadronization, which can no 
longer be treated perturbatively, since the hadrons are formed by the purely nonper- 
turbative effects of color confinement. Therefore the dynamics of confinement should 
be considered in the transport descriptions. 

There exist different approaches to deriving a transport theory that is based on an ef- 
fective quark model incorporating the dynamics of the nonperturbative QCD vacuum. 
Zhang and Wilets P| have derived a transport theory based on the Nambu Jona- 
Lasinio model in order to estimate chiral symmetry effects in heavy-ion collisions; 
this is conceptually close to our work. However, besides the absence of confinement, 
these authors do not actually perform dynamical simulations. Also soliton-soliton 
collisions in the framework of the Skyrme- model and of a Skyrme-like a-model 

have been performed. However, these calculations do not include any explicit 
quark degrees of freedom. 

We assume that the nonperturbative effects can be modeled by the exchange of a 
scalar field, as it is successfully done in the Friedberg-Lee model (for a review see 
Top. Since in this model the hadron surface is generated dynamically, it is well 



suited for dynamical simulations. 

In recent articles |Tl], |12[ we presented the derivation of the transport theory based on 
the Friedberg-Lee model and applied it to simulations of nucleon-nucleon collisions 
for bombarding energies of the order of a few GeV. In this work color degrees of free- 
dom were neglected. Here we present the derivation and first dynamical applications 
of a transport model with color confinement. After having fixed the strong coupling 
constant with the help of the string constant, we will investigate an excited state of 
the nucleon, the color isovector oscillation mode, showing the dynamical impact of 
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color confinement on the quark motion. 
The phenomenological Lagrangian of the Friedberg-Lee model reads [|T3|, P 



L = ^{t^,d^ - mo - goa)^ + ^{d^crY - U{a) - ^K(a)F;,Fr - ig.^i,^^A>:^ , (1) 



where \E' denotes the quark fields, a is the color singlet scalar field representing 
the long range and nonabelian effects (of multi gluon exchange), and the last term 
contains the interaction of the residual abelian (one gluon exchange) color fields 

F^u = - duAl (2) 

All the nonabelian effects are assumed to be absorbed in the color dielectric function 
K{a) which is chosen such that k vanishes as a approaches its vacuum value a = 
outside the bag and k = 1 inside. U{a) is the self-interaction potential for the scalar 
a-field containing cubic and quartic terms. The Euler-Lagrange equations derived 
from dll) are 

{^^^{id, - ig.^AD - mo - ^?o^)^ = , (3) 
d.d^a + U'{a) + i/s:'(a)F;,Fr + Qo^^ = , (4) 

5'^W=J:, (5) 



where (^ is the Dirac equation describing the motion of the quarks, (^) is the Klein- 
Gordon equation for the time evolution of the cr-field, and (^ are the Maxwell equa- 



tions for the gluon fields with the color charge currents = —ig^'^'yy^^^!. Solutions 
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of these coupled set of nonlinear equations have up to now only been obtained in 
static calculations |T^, since dynamical simulations suffer from instabilities, 

due to the small current quark mass and large Dirac-sea contributions. Here we will 
determine the dynamics of the quarks by a transport equation derived by Elze and 
Heinz [IITI . These authors have shown on a semiclassical basis, that the following set 



of equations for the phase space distribution functions in the abelian limit can be 
derived: 

(p^9^ - m*d,m*d>;)f\x,v) = gvP^.FrdPnx,p) , (6) 



{p,d^ - m*d,m*d^)nx,p) = g.p,Frdtf{x,p) . (7) 

Note that these equations (§) and (|^ are the transport equations of motion from 
the Friedberg-Lee Lagrangian, whereas the corresponding equations (2.14) of [|T^ are 
derived from the QCD Lagrangian. /° and 

f{x,p)= j dQf{x,p,Q) , r{x,p) = JdQQ^f{x,p,Q) (8) 

are the color singlet and the color octet distribution functions of the general distri- 
bution function f{x,p,Q) for colored particles, with the space-time vector x, the 
four-momentum p and an eight component c-number color vector, Q = Q"^. From 
these distribution functions we obtain the scalar density ps = ^"^y and the color 
charge density Jq needed for the equations of motion (|) and (j^) by 

with the effective mass m* = mo + g^a and the factor 77 = 4 from spin and flavor 



degeneracy |[T2] 



Because the color fields in (Q) and (^) are purely abelian the equations decouple 
to eight identical Maxwell equations. Adopting the Coulomb gauge V{kA) = 0, these 



equations can be written as 18 



V{kVAo) = -Jo . (10) 



Following the general assumptions of bag models we assume that the contribu- 
tions from the current is already absorbed in the current quark mass mo- Hence 
we focus on the colorelectric field described by (|TUp only, especially in view of qq- 
production from strong colorelectric fields. 
We make the following particular ansatz for and 

fix,p) = fix,p) + fix,p), (11) 
f\x,p) = ... = f{x,p)=0, (12) 
f{x,p) = fix,p)-fix,p). (13) 



This ansatz corresponds to working in a U(l) subspace of SU(3)-color. This can be 
seen by remembering that there always exists an (abehan) gauge transformation such 
that the color vectors of the quarks contribute only to the commutating generators 
of SU(3), and [^. Adopting this expansion, the quark (/) and anti-quark (/) 
distribution functions appearing in (|Tl|) and (|T3]) can be interpreted as quark and 
di-quark distribution functions. Nevertheless, we will in the following refer to (/) 
as the anti-quark distribution function, the U(l) equivalent of the SU(3) di-quark. 
Inserting the ansatz (|T1|)-(|TB|) into the equations for the color moments (^) and (|^ 
we end up with 

ip.d^- - m*id,m*)d;)fix,p) = g,p,F^''dlf{x,p) (14) 
{p^d^^ - m*{d,m*)d^)f{x,p) = -g.p.F^^'dlJix.p) , (15) 

which is a set of usual Vlasov equations describing the motion of charged particles in 
a selfconsistently generated scalar and vector field, determined by equation (H) and 
(p!0|), respectively. The coupling of quark and antiquark degrees of freedom is there- 
fore only provided by their interaction with the mean fields a and the colorelectric 
field. 

Let us first examine the static limit of the transport equations. From quantum me- 
chanics we know that the color charge density in a hadron has to vanish locally 



21 



I' 



dMf{x,p) - f{x,p)) = f,{x) =< N\Q\x)\N >= , (16) 

because the nucleon ground state | > is a color singlet and a color octet operator. 
From that we conclude that / = / and hence there is no colorelectric field in the 
groundstate. As shown by Vetter et. al |T^, the distribution functions have to be of 



a local Thomas Fermi type 

f{x,p) = f{x,p) = eifi-Lu), (17) 



with the Fermi energy /i. After inserting the local Thomas Fermi distributions (|Tj 
into the integral expressions (^ the fermions can be integrated out, so that we are 



left with determining the sohton solution of the remaining cx-field equation (^. This 
solution is then obtained by standard methods used by Vetter et. al [|T^. Comparing 



this solution with grounstate properties of the nucleon, the parameters of the model, 
go and the ones of the potential U{cr), can be fixed. 

Since the colorelectric field vanishes in the static (ground state) solution (0), the 
strong coupling constant a, = gl/4:7c cannot be fixed from ground state properties. 
In order to determine its value we investigate the colorelectric energy in the cavity 
by forming a string like configuration. The string tension r, i. e. the coefficient of 
the linear rising gg-potential 

k , 4ao 



Vgg{r) = h rr , k 

r 



within soliton bag models has already been determined by several groups |^5|, |I9| ^ 
In these studies, the quark motion is neglected and an instantaneous electric field 
between them is assumed ab initio. We proceed differently and calculate the response 
of the colorelectric field to an adiabatic separation of the quark and anti-quark dis- 
tributions with a constant velocity. Since we work in the testparticle picture p3|, a 
fixed initial velocity Vz is assigned to each particle belonging to / and —Vz to each 
particle belonging to /. The velocity of the quarks is kept fixed during the time 
evolution. In this way the source densities for the time evolution of the cx-field and 
the instantaneous A°-field equations (|) and (|T0|) are obtained. Several values of Vz, 
ranging from Vz = 0.1c to Vz = 0.5c, have been applied to assure the independence 
of the colorelectric field on Vz- 

At the beginning of the motion we see in fig. [l| how the cavity is deformed by the 
motion of the quarks {t = 3fm). While the charge distributions move apart the 
colorelectric field, as a solution of (p!OD, builds up (fig. |^). The genuine feature 
of this field is, that it is parallel everywhere at the boundary of the cavity, which 
is due to the von Neumann boundary conditions, forcing the normal component of 
the displacement field D = kE to vanish. Hence we find the typical string like be- 
haviour of the i?-field in the later stages of the time evolution (fig. Qb), namely being 



constant and almost independent of the axial coordinate in the region between the 
quarks. Determining the colorelectric energy Egiue = | / d?rDE as a function of the 
gg-distance provides us with a measure for the string constant r. Depending on the 
value of Us the string constant is shown in fig. 0. In order to obtain the generally 
accepted value of r ^ IGeV / fm, we need to have values of as ~ 2, which is in 
agreement with the values of MIT-bag and other Friedberg-Lee model calculations 
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Having fixed the strong coupling constant, we now perform the time evolution of 



the full model equations for the a-field (P, the gluon field (10) and the transport 



equations (|T^ - (0). As a first example, we excite the soliton ground state at time 
t = by shifting the equilibrium (Thomas Fermi) momentum distribution (|T^ of the 
quarks with a fixed momentum +pz and that of the anti-quarks with —pz- (Values 
used for pz range from pz = fi = 350MeV to pz = 30/i). A self-consistent time 
evolution of the equations of motion (^ and (|T3p-(^) is then generated by a numer- 
ical integration sceme as used by Vetter |jl2|. The instantaneous colorelectric field is 
determined in every timestep of the evolution by a finite element method developed 
by Mitchell . The general features of the resulting process can be summarized as 
follows: 

Depending on their initial momenta, the quarks move apart along the z-axis and the 
electric field builds up. As a consequence, the quarks lose their kinetic energy to 
the E-fie\d and additionally to the deformation of the cavity to a string like con- 
figuration. When the entire kinetic energy is deposited into the E- and a-field, the 
motion is stopped and starts to reverse. The quarks then pass through each other 
at the origin and the process is repeated in the opposite direction, leading to an 
oscillation of the gg-system. If we follow the process further, the coherent motion of 
the gg-system is disturbed due to an increasing chaotic motion of the testparticles 
inside the deformed cavity As a consequence, the strength of the colorelectric 
field decreases as shown in (fig. ^) for an initial energy of 3GeV, accompanied by a 
steady elongation of the cavity. This behaviour continues until the total loss of the 
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colorelectric energy is used for a large amplitude vibration of the bag (fig. For 
very large times the superposition of the q and q distributions is complete and the 
colorelectric field vanishes, so that the final state is again a locally color neutral one 
as it was at t = 0. 

Note that the elongation of the cavity shown in fig. |^ is caused by the initial mo- 
mentum pz that acts as an additional pressure on the cavity. A further effect of the 
chaotic motion to be expected is the damping of multipole modes in Vlasov-type 
transport models leading to an equal distribution of the particle momenta to 
all degrees of freedom. This damping can hardly be seen in our particular scenarios, 
since the excitation energies (typically of the order of a few GeV) are too large, so 
that the equipartition of the particle momenta happens on a much larger time scale. 
For lower exitation energies of ~ SOOMeV, however, this one body friction effect, 
driving the bag to a spherical shape, can be seen. 

Summarizing the statements above, we can say that independent of the relative mo- 
mentum of the quarks, the motion is always restricted to the cavity by the mutual 
interaction of the ii^-field and the cx-field with the quarks inside the dynamically 
deformed cavity. Thus, the inherent feature of the Friedberg-Lee model, color con- 
finement, is maintained in our dynamical model, i. e. for an arbitrary momentum of 
the quarks. 

Further applications of our model will be the effects of color gluon fields in nucleon- 
nucleon collisions, similar to our investigations in the colorless model |Tl|, |12|, or the 
description of quark gluon plasma scenarios, in particular the dynamical treatment 
of the hadronization process. Also a consistent description of gg-production can be 
implemented via the Schwinger mechanism. Work along this direction is in progress. 
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Figure 1: Baryon density p — /(/ + f)d^p (left) and the u-field (right) at time 

t = 3/m (upper) and t = 40/m (lower). The equidistant lines of constant value start 
at 0.1/m~^ for the quark density p with an increase of 0.1/m~^. The values for a 
start at — 0.125/m~^ with an increase of 0.05/m~^. 



Figure 2: The electric field E inside the cavity, a) At time t — 3/m/c and b) at time 
t — 40/m/c. Note, that b) has a rescaled z-axis to fit the page. 



Figure 3: String constant r as a function of the strong coupling constant Ug- 



Figure 4: Time evolution of the colorelectric energy up to i = 100/m/c resulting from 
an exitation energy of 3GeV (dotted line). For the solid line the strong oscillations 
have been averaged out. 



Figure 5: Expectation value of the 2;-coordinate < z >— J J zp{x,p) dPx d^p as a 
function of the time t resulting from an exitation energy of 3GeV (dotted hue). For 
the solid hue the oscillations have been averaged out. 
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